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We use a 15 lattice QCD heavy-light wavefunction to generate a single 
parameter, model independent description of the Isgur-Wise function. Using 
recent data we find the zero-recoil slope to be £'(1) = —1.16 ±0.17, while the 
second derivative turns out to be £"(1) = 2.64 ± 0.74. 
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1 Introduction 



The recent development of the Heavy Quark Effective Theory [||] yields an expression 
for the B — > D^Wi decay rate in terms of a single unknown form factor, the Isgur- 
Wise function (IW). This function is absolutely normalized at zero recoil point 
up to corrections of order 0. It is currently believed that these corrections 
can be calculated with less than 5% uncertainty ||, which would allow a precise 
determination of the CKM matrix element \V c b\ from the study of B — > D^Wi 
decay SIS db function of the D^*> meson recoil. 

In this note we take a somewhat different approach. We use an established 
LQCD result |4], |5| for the heavy-light wavefunction to express the IW function 
in terms of a single parameter, and then constrain that parameter using recent 
improved experimental results by the CLEO II group || in exclusive semileptonic 
B -> D^Wi decay 0. 

2 Numerical procedure and results 

If one knows the heavy-light meson wavefunction in its rest frame, and the energy 
of light degrees of freedom E q , then the IW function for a given u = v ■ v' (where v 
and v' are 4- velocities of two hadrons), can be computed || from 

««) = ;^r<W^T r)> ■ (1) 

where 

(A) = / dr r 2 R(r)A(r)R(r) . (2) 
Jo 

From the above the first and second derivatives at the zero recoil point are [§, [| , 

= -(^ + ^V>), (3) 

= l + l E "( r2 ) + YE E t( r4 ) ■ ( 4 ) 
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Higher derivatives could similarly be computed if desired. 

The heavy-light wavefunction has been recently computed by a lattice simulation 
In order to use the above formulas and the LQCD "data", we need an analytic 
expression for the wavefunction. Instead of trying to find some specific analytic form 
that would describe the behavior of the lattice data close to the origin and at large 
r, we expand the lattice wavefunction in terms of a complete set of basis functions. 
We then truncate the expansion to the first N basis states, hoping that we are be 
able to find a good description with a small number of basis states. In other words, 

N-l 

R ls (r) ~ Yl ^ ■ ( 5 ) 

i=0 



The quasi-Coulombic basis set |T(J which we have chosen is particularly well suited 
for relativistic hadron systems, and it is given (for s-waves) by 



Ri(r) 



e-t»Li(2Pr) , (6) 



\(i + 2)(i + l) 

where L\ are associated Laguerre polynomials and (5 is a scale parameter. Substi- 
tuting (||) and (|6|) in the expression ([[J) yields 

p( \ 2 (i + 2)!(j + 2)! 

u + 1 i=o 3 -=o ^/(i + 2)(i + l)(j + 2)(j + l) 



where 



^ _ ^ ^ (-l) m+n (m + n + 1)! sin [(m + ra + 2) arctan (a)] ^ 
m =on=o m\n\(i - m)\(m + 2)!(j - n)\(n + 2)!a(l + a) m a" 



and 



(9) 



/? V w + 1 ' 

In deriving this expression we have used the Laguerre polynomial representation 

L?(x) = ± ( ' + Q) ; . (10) 

^T m! (z — m)!(m + a)! 
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Our radial basis states depend on the meson size parameter /3, which we estimate 
from the exponential falloff of the lattice data (~ e _/3r ) to be somewhere between 
0.35 and 0.45 GeV. Once (3 has been fixed, we vary the coefficients q from (^) in 
order to best fit the lattice data. Fortunately, for any j3 within the expected range 
we are able to find an excellent approximation to the lattice wavefunction with only 
3 basis states, with \ 2 °f about 1.1 per degree of freedom (we have assumed an 
uncertainty of 10% in the wavefunction values JOJ). Since all fits were essentially 
equivalent, and the wavefunctions were nearly identical up to r = 20 GeV' 1 , we 
have chosen the intermediate value (3 = 0.40 GeV. We emphasize though that none 
of the details of the basis function representation are important to our final result. 
In Figure [I] we show our radial wavefunction for the IS state, together with LQCD 
data points, for (3 = 0.40 GeV. As one can see from the figure the agreement 
is excellent, even though there is some ambiguity in the data, especially at large 
r, where the effects of the small lattice size become large. The three basis state 
coefficients (ft = 0.40 GeV) are found to be 

c = 0.9985 , 

ci = -0.0221 , (11) 
c 2 = 0.0500 . 

Unfortunately, the energy eigenvalue associated with the lattice wavefunction is 
not easily interpreted. The value of E q depends sensitively upon the lattice spacing 
H [TT]] . Therefore, in the calculation of the IW function we treat E q as a parameter. 

From (0) and (|TTD we can now compute in terms of E q . In Fig. ^|we show 



the IW function for the range of 

E q = 0.306 ± 0.040 GeV , (12) 

which corresponds to a one standard deviation corridor for the seven CLEO II || 
data points (or \ 2 of about 0.65 per degree of freedom). We also show (full line) 
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the IW function corresponding to the best fit. Finally, we use (^) and (|4]) with the 
allowed range of E q ( [T2"| ) to evaluate the first and second derivatives at zero recoil 
point. The results are 



2.64 ±0.74 . 



-1.16 ±0.17 , 



(13) 
(14) 



The best fit corresponding to E q = 0.307 GeV, £'(1) = -1.15, and £"(1) = 2.56, has 
X 2 of 0.38 per degree of freedom. 

3 Conclusion 

Since the IW function is non-perturbative any parametrization necessarily must 
contain some physical input. We have used a reliable lattice result for the heavy- 
light meson wavefunction to compute the IW function in terms of a single parameter 
E q . By comparing with experimental data we determine the allowed range of this 
parameter. Among the direct uses of this parametrization is a more believable 
extraction of the IW slope at zero recoil point and for the first time a reliable 
estimate of the second derivative. Previous slope estimates have assumed £(w) can 
accurately approximated by £(co) ~ 1 ± £'(l)(co> — 1), but this inevitably leads to an 
overestimated (not sufficiently negative) value for the slope 0. 
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FIGURES 



Figure 1: S-wave radial wavefunction obtained from the fits to the LQCD data. We 
used (3 = 0.40 GeV and N = 3 basis states. The wave function is normalized as in 



Figure 2: Two limiting cases of the IW function calculated from the IS ra- 
dial wavefunction shown on Fig. |], corresponding to E q = 0.346 GeV and 
E q = 0.266 GeV^dashed lines). The best fit to the data is for E q = 0.307 GeV 
(full line). For the sake of clarity, the error bars are shown only for the CLEO II 
data. 
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Figure 1 
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Figure 2 



9 



